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In this talk we will discuss the physics reach of the atmospheric neutrino data collected by a future 
Qh megaton-class neutrino detector. First we will discuss the potentialities of atmospheric neutrinos 

Q_i on general basis, presenting our results in the form of neutrino oscillograms of the Earth: contour 

plots in the neutrino energy-nadir angle plane. In this context we will analyze in detail the various 
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signatures related to 613, to the neutrino mass hierarchy, to the octant of 023 an d to the 8 C p phase 
which appear for different values of the neutrino energy and baseline. Then we will consider 
concrete experimental setups, showing that synergic effects exist between atmospheric and long- 
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baseline neutrino data: the combination of the two data sets is much more powerful than the 
simple sum of the sensitivity of each individual data sample. 
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1. Introduction 

Despite their pioneering contribution to the discovery of neutrino oscillations, it is in general 
assumed that atmospheric neutrino data will no longer play an active role in neutrino physics in 
the coming years. This is mainly due to the large theoretical uncertainties arising from the poor 
knowledge of the atmospheric neutrino fluxes, which strongly contrast with the requirement of 
"precision" needed to further enhance our knowledge of the neutrino mass matrix. In this talk, 
we will show that despite these large uncertainties the atmospheric neutrino data collected by a 
megaton-class detector will still provide very useful information. 

The main strength of atmospheric data is its very broad interval in neutrino energy (E v ) and 
baseline (determined by the nadir angle V ). In order to provide a global view of this whole range, 
in this section we will make extensive use of neutrino oscillograms of the Earth, i.e. contours plots 
in the neutrino energy-nadir angle plane [1,2]. Consider a bin centered at (0 V , E v ) with extensions 
A0 V and AE V . The numbers of expected (NnJ and observed (./Vex) events in such bin are of course 
proportional to its size. Hence its contribution Ax 2 to the total x 2 function is: 



As we will see, different regions of the (E v , V ) plane will exhibit characteristic structures whose 
position and size is determined by various neutrino parameters. The function t; provides an easy 
way to highlight which of these regions mostly contribute to the x 2 - Hence in the following we 
will present isocontours of t; . 

2. Sensitivity to oscillation parameters 

It is well known that for non-zero value of $13 matter effects induce a resonance in the — > v e 
conversion probability at E v ~ 3 4- 6 GeV. The precise position of this peak in the (E v , V ) plane 
depends on the value of $13, so that in principle atmospheric neutrinos could be used to measure this 
angle as long as it is larger than about 3°. However, in practice the sensitivity is limited since the 
Vfi — ► v e signal is unavoidably diluted by the v e — > V e background, and also since the atmospheric 
neutrino flux at E v ~ 6 GeV is considerably suppressed. Therefore, although some sensitivity is to 
be expected in a megaton detector, it is likely that atmospheric neutrinos will not be competitive 
with dedicated long-baseline and reactor experiments in the determination of 613. 

The sensitivity to the hierarchy is illustrated in Fig. 1. We can distinguish two effects: 
- in the v e channel there is a peak arising from the high-energy resonance just discussed for 613. 
This signature is only visible if 613 is large enough. Note that in order to determine the hierarchy 
it is not sufficient to see the resonance: it is also necessary to tell whether it occurred for neu- 
trinos (normal hierarchy) or antineutrinos (inverted hierarchy). It is therefore crucial to have a 
detector capable of charge discrimination. In the case of a non-magnetic detector such as a Wa- 
ter Cerenkov, a statistical separation of v and V events can still be obtained by measuring both 
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Inverted Hierarchy 
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Figure 1: Sensitivity to the neutrino mass hierarchy. Left: sin 29^ — 0.125 (N ex ) versus 613 = (N t h). 
Right: normal hierarchy (A^x) versus inverted hierarchy (N t h). The undisplayed parameters are Ato 2 j = 



8x 1(T 5 eV 



, Am^j =2.5 x 10" 



eV 2 , tan 2 Q n = 0.45, sin 2 2 3 = 1 and <5 CP = 0. 



single -ring and multi-ring events, since the relative contribution of neutrinos and antineutrinos 
to these two data sets is different. 

- in the channel we observe a characteristic structure of peaks and deeps, which in principle 
can be used to obtain information on the mass hierarchy. However, due to the fast alternation 
of positive and negative regions a very good resolution on both neutrino energy and direction is 
needed, otherwise the signal averages to zero and all information is lost. Furthermore, the posi- 
tion of the structures is affected also by other oscillation parameters, so that a poor knowledge of 
these parameters will result in a considerable loss of sensitivity. Note that this signature is also 
present for 0^ = 0, however it is strongly suppressed and hardly detectable. 

In brief, atmospheric neutrino data can provide useful information on the neutrino mass hierarchy, 
but only if 613 is sufficiently large and a detector with good resolution and charge discrimination 
capability is available. 

The sensitivity to the octant is one of the topics where atmospheric neutrino are mostly useful. 
As can be seen in Fig. 2, we have two characteristic signatures: 

- at low energy (E v < 1 GeV) we observe an excess of e-like events and a deficit of /x-like events 
for 023 < 45° with respect to 023 > 45°. The v e signal is about four times larger than the 
one [3]. This effect is due to subleading oscillations induced by Am^j, and is present also for 
013 =0. For 0n 7^ the neutrino flux arriving at the detector is modulated with the very fast 
Ara^ oscillations, but the effect persists on average. Since this signature appears with the same 
sign for neutrinos and antineutrinos, no charge discrimination is required for its identification. 

- at high energy (E v > 3 GeV) we observe a deficit of e-like events and an excess of /i-like events 
for 023 < 45° with respect to 023 > 45°. This effect is again related to the matter resonance 
discussed for 0i3 and the hierarchy, and indeed it appears only for ©13 / and only for v's or 
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Figure 2: Sensitivity to the octant, for 0n = (left) and sin 2 ©13 = 0. 125 (right). Both plots show sin 023 = 
0.4 (A^x) versus sin 2 023 = 0.6 (A^h)- The undisplayed parameters are Am^ = 8 x 10~ 5 eV 2 , Am^ = 2.5 x 
1(T 3 eV 2 , tan 2 0i 2 = 0.45 and <5 CP = 0. 



v's. Therefore, in a detector without charge discrimination this signal will be suppressed. 
Note that the presence of a low-energy effect independent of $13 guarantees a minimum sensitivity 
to the octant from atmospheric neutrinos, provided that the deviation of 623 from maximal mixing 
is large enough. This is a unique feature which will prove very synergic with long-baseline data, as 
we will show in the next section. 



3. Synergies with long-baseline experiments 

So far we have discussed the potentialities of atmospheric neutrinos in general terms. Let us 
now consider three concrete experimental setups [4] : 

- a Beta Beam (j8B) from CERN to Frejus (130 Km). We assume 5 years of v e from 18 Ne and 5 
years of v e from 6 He at 7 = 100, with an average energy (E v ) = 400 MeV. For the detector we 
assume the MEMPHYS Water-Cerenkov proposal, corresponding to 3 tanks of 145 Kton each; 

- a Super Beam (SPL) from CERN to Frejus (130 Km). We assume 2 years of Vu and 8 years of 
Vn running, with an average energy (E v ) = 300 MeV. Again we use MEMPHYS as detector; 

- the T2K phase II (T2HK) experiment, corresponding to a 4MW super beam from Tokai to 
Kamioka (295 Km), with 2 years of and 8 years of v^. The detector is the proposed Hyper- 
Kamiokande, rescaled to 440 Kton for a fair comparison with the /3B and the SPL. 

As widely discussed at this conference, a characteristic feature in the analysis of future LBL 
experiments is the presence of parameter degeneracies, which pose a serious limitation to the deter- 
mination of 0i3, 8 C p and the sign of Ara^. This problem is illustrated in Fig. 3, where we show the 
allowed regions in the plane of sin 2 2 613 and 5 CP for /3B, SPL and T2HK. As visible in this figure, 
for the /3B the intrinsic degeneracy cannot be resolved, while for the super beam experiments SPL 
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95% CL regions for the (hW), (H tr O wr ), (H wr O tr ), (H wr ( solutions 
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Figure 3: Allowed regions in sin 2 20b and 5cp for LBL data alone (contour lines) and LBL+ATM data 
combined (colored regions) [4]. H tr/wr and O tl/wl refers to solutions with the true/wrong mass hierarchy 
and octant, respectively. The true parameter values are AotIj = 7.9 x 1(T 5 eV 2 , Am^ = 2.4 x 1(T 3 eV 2 , 
sin 2 0i2 = 0.3, sin 2 20B = 0.03, sin 2 23 = 0.6 and <5 CP = -0.857T. 

and T2HK there is only a four-fold degeneracy related to the sign of Am 2 j and the octant of 023- 
Once atmospheric data are included in the fit all the degeneracies are nearly completely resolved, 
and the true solution is identified at 95% CL. This clearly show the presence of a synergy between 
atmospheric and long-baseline data: at least for this specific example, the combination of the two 
sets is much more powerful than the simple sum of each individual data sample. 

4. Conclusions 

In this talk we have discussed the potentialities of atmospheric neutrino data in the context 
of future neutrino experiments. We have shown that despite the large uncertainties in the neutrino 
fluxes atmospheric data will still provide useful information on the neutrino parameters, due to 
their very broad range in neutrino energy and nadir angle. In particular, we have proved that the 
sensitivity obtained by a combination of atmospheric and long-baseline data is much stronger than 
the one achievable by each data set separately. 

Work supported by MICINN through the Ramon y Cajal program and through the national 
project FPA2006-01105, and by the Comunidad Autonoma de Madrid through the HEPHACOS 
project P-ESP-00346. 
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